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Abstract 



Quantum tunneling between two potential wells in a magnetic field can be 
strongly increased when the potential barrier varies in the direction perpendic- 
ular to the line connecting the two wells and remains constant along this line. 
A periodic structure of the wave function is formed in the direction joining 
the wells. The resulting motion can be coherent like motion in a conventional 
narrow band periodic structure. A particle penetrates the barrier over a long 
distance which strongly contrasts to WKB-like tunneling. The whole problem 
is stationary. The coherent process can be influenced by dissipation. 
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Can a particle move under a long and almost classical potential barrier where classical 
motion is impossible? 

According to Wentzel, Kramers, and Brillouin (WKB) [1], there is a finite probability 
w ~ exp(— A) of quantum tunneling through an one-dimensional potential barrier. This 
probability becomes negligible for semiclassical barriers when A = 2\mS/ h and the classical 
under-barrier action ImS* is big. In two- dimensions the most convenient way to calculate 
the exponent A is by use of a classical trajectory x(t), y{r) in imaginary time t = it [2-6]. 
The trajectory goes in a classically forbidden area (under the barrier) and connects two 
classically allowed regions having zero velocities at the borders of classical regions. The 
classical action, constructed by means of this trajectory, is called Euclidean action and 
determines WKB-type exponent A. The method of classical trajectories in imaginary time 
is powerful and relatively simple since it allows to determine a tunneling probability in 
the main (no pre-exponent) approximation exp(— A) just only solving Newton's equation of 
motion. 

The problem of quantum tunneling in magnetic field was addressed in Refs. [7-9] . In the 
Landau gauge there is a parabolic gauge potential mu 2 (x — xo) 2 /2 superimposed upon the 
tunnel barrier potential. The cyclotron frequency is uj c = eH/mc and tunneling occurs in 
the x-direction. If the tunnel barrier is not a constant, containing weak impurity centers, 
x Q becomes spatially dependent resulting in a variable gauge potential of a sawtooth shape 
instead of pure parabolic one. This potential is "pinned" by impurities, separated by the 
characteristic distance b, and repeats their positions [7,8]. In the regime of a strong magnetic 
field, when the energy muj 2 b 2 /2 exceeds a height of the tunnel barrier, an electron tunnels 
incoherently through each peak of the gauge potential. For this incoherent motion the total 
probability of tunneling is a product of partial ones. 

This process can be elegantly described in terms of classical trajectory in imaginary 
time (t = it) [10] (see also [11,12]). The classical trajectory is constructed in the following 
way: x(t) is real and corresponds to a translational motion in the direction of tunneling 
(a;- axis) but the transverse coordinate y = it]{t) is imaginary and performs oscillations in 
t. The oscillations are produced by an oscillatory force in the x-direction since impurities 
are distributed along the x-axis. This potential force is balanced by the ^-component of the 
Lorentz force muj c drj/dr. The tunnel potential U(x,y) is supposed to be even with respect 
to y and the classical trajectory satisfies the relation U [x(t), ii](T)] = E. At the ends of the 
trajectory the relation 77 = holds which allows to match with the physical y = 0. 

If the magnetic field is not big, kinetic energy enters the game and a scenario of tunneling 
can be dramatically different. We consider the two-dimensional tunnel potential which in 
the barrier region has the form U(x,y) = u(y), tunneling occurs in the x-direction, and the 
magnetic field is aligned along the z-axis. In this case a variable gauge potential can form 
under the barrier a certain structure periodic in the direction of tunneling with a period 
Ax. An electron moves in this periodic potential in the way similar to conventional motion 
in a periodic structure with a narrow energy band AE where tunneling processes through 
subsequent periodic barriers are strongly coherent. For a conventional periodic structure a 
wave packet can pass over a long distance with no exponential decrease in amplitude but 
only with a delay time h/ AE as a velocity in a narrow band is proportional to AE [13]. In 
our case, since the potential is ^-independent in its barrier part, a period in the x-direction 
Ax = y/2\E \/m At (^2\E\/m plays a role of velocity) has an intrinsic nature. As shown 
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below, the time At is a period of oscillations in the well associated with the potential 



v{rj) = u{irj) 



muoz 



(1) 



The particle energy E is negative, rj = ^2\E\/muj%, and the minimum of u(y) corresponds 
to u(0) = 0. The transverse oscillatory motion in the direction of y = irj is coupled to the 
translational motion (x-direction) due to the Lorentz force. 

When the magnetic field is close to the certain value Hr the probability of tunneling 
w ~ exp(— A) through a long barrier becomes not exponentially small, like for a conventional 
narrow band dynamics, and corresponds to A — > 0. This is a situation of Euclidean resonance 
studied in papers [14-18] for tunneling through nonstationary barriers when also A — > at 
a certain value of an ac amplitude. Therefore, a phenomenon of Euclidean resonance has 
rather a general nature since it occurs also in a static barrier. As argued in this paper, 
the magnetic field sets a long distance under-barrier coherence which allows under-barrier 
motion over a long distance. This strongly contrasts to WKB-like tunneling. 

Below we formulate the above arguments in terms of a classical trajectory in imaginary 
time (t — it). A potential barrier, part of which is plotted in Fig. 1, is even with respect 
to y, it does not depend on x at < x < R where U (x,y) = u(y) (u(0) = 0), and the 
function U (x, 0) has jumps at x = and x = R. Tunneling occurs between two classically 
allowed regions x < and R < x. For convenience, the potential in Fig. 1 is drawn in a 
way that it is a constant at x < and R < x. This condition is not necessary and U (x, y) 
can correspond to tunneling between two quantum wires or two quantum dots. In order 
to calculate a tunneling probability in the exponential approximation one can know only a 
classical trajectory in imaginary time connecting two points {x = 0,y = 0} at the moment 
t = t and {x = R, y = 0} at the moment r = 0. Classical equation of motion under the 
barrier have the form 



d 2 x drj 
m —r + mu c — = 0; 
or 1 or 

The total energy conserves 



d 2 i] dx duiirj) 

m + mw c — H — 

or 1 or or] 







(2) 



E 



The conditions to Eqs. (2) are 

drj 

According to Eq. (3), 



(IK 



dr] 



TO 



(3) 



= 0; 77(0) = 77(7-0) = 



(4) 



dx 



dx 



TO 



2\E 



m 



(5) 



For a smooth potential all velocities should be zero at the ends of a trajectory but in our case 
dx/dr is finite due to jumps in the potential energy at x = and x = R. The differential 
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equations (2), with respect to the functions dx/dr and rj, depend on three parameters which 
should be determined from the six conditions (4) and (5). Since the functions dx/dr and r\ 
are periodic among the six conditions there are only three independent ones. 

The probability of tunneling w ~ exp(— A) from x = to x = R in Fig. 1 is expressed 
through the Euclidean action [2,3,5,6,10-12,19-21] 



TO 



A=% I dr 



m ( dx\ 2 m ( di]\ z dx 



2\aF) -^Wr) + ™ c r ] -+u(ir ] )-E 



(6) 



Without a magnetic field rj = and the action (6) coincides with the WKB expression. 
With the solution of the first equation (2), dx/dr = —uo c (rj + 770), the action (6) reads 
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A = - / dr 
ft Jo 



TO 



m ( di]\ 2 dx 



(7) 



By means of Eq. (3) the expression (7) takes the form 
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TO 



A = A WKB --J dr[E-v (77)] (8) 

where Awkb = 2 v /2m\E\R/h comes from the last term in Eq. (7) and it is the WKB action 
related to one-dimensional tunneling across a rectangular barrier of the length R with the 
energy \E\ below the barrier top. Awkb in Eq. (8) is generic with the conventional under- 
barrier action in a multi-dimensional case [2,3,5,6]. The second (negative) term in Eq. (8) is 
solely due to the magnetic field and corresponds to the transverse motion. The total energy 
now has the form 

In contrast to the x-motion, kinetic energy in the transverse channel does not change sign 
compared to the physical trajectory (dy/dt) 2 = (drj/dr) 2 . Therefore, the transverse motion 
occurs at the region where v (77) < E and the second term in Eq. (8) is negative. With the 
expression (9) one can describe an oscillatory motion in the potential v{rj). We specify a 
shape of v{rj) as in Fig. 2 where v (At/) = E. The trajectory rj(r) is drawn in Fig. 3(a) where 
the period Ar, according to Eq. (9), is 

Ar = V2^l - dV (10) 

Jo ^E - v( V ) 

The trajectory x(r) is shown in Fig. 3(b). Each cycle of r/(r) in the potential well in Fig. 2 
results in the translation of x{r) by Ax determined by 

^= (11) 
VE ~ v( V ) 

The trajectory in Fig. 3 looks qualitatively similar as one related to tunneling through a 
barrier slightly violated by impurities in a strong magnetic field [10,11]. The difference is 
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that the oscillations in Fig. 3 have an intrinsic nature but in [10,11] they are determined by 
distributed impurities. 

This method of trajectories is applicable when the distance between the wells is R = NAx 
(and also r = TV At), since at the ends (x — and x = R) the condition rj = should hold. 
Fig. 3 is plotted for N — 3. By means of (9), the last term in Eq. (8) can be written as 
(-NAA) where 

AA = -I— / driy/E-virj) (12) 
Jo 

and the action (8) takes the form 

At a small magnetic field E — v(irj) ~ —u{irj) and Ax ~ ln(l/iJ) is logarithmically big. 
Therefore, at H — > the action (13) turns to its conventional limit A WKB . 

The action (13) and its WKB part are shown in Fig. 4(a) at the points R = NAx where 
they have only sense. If R is not NAx another method of calculation is required. The action 
A is substantially reduced compared to its WKB part A-^kb- When the magnetic field is 
close to the certain value Hr the action (13) is A ~ (Hr — H) and may turn to zero. This 
means that the tunneling probability w(H) becomes not exponentially small at H = H R . 
As mentioned above, the phenomenon, when A — > 0, is called Euclidean resonance. 

The method of trajectory used corresponds to one-instanton approach and it holds at 
H < Hr when exp(— A) is small. At H > H R the exponent exp(— A) is not small and one 
should apply a multi-instanton approach which accounts all powers of the exponent. This 
is a matter of a further study. One has to expect that at H > H R the tunneling probability 
w(H) also decays in a manner as at H < Hr having the peak at H = Hr plotted in Fig. 5. 
The peak width, as one can show, is roughly Hr/Awkb- The condition R = NAx(H), 
when the above method of trajectories is applicable, holds at certain values of the magnetic 
field H = h-N shown in Fig. 5. The relation 

w(h N -i) = w(h N )exp(-AA) (14) 

says that always exists some for which w^h^) is not smaller than exp(— AA). 

The oscillatory structure of the trajectory in the x-direction indicates analogous oscil- 
lations in the wave function. At H < Hr and R = NAx the function \ip(x, 0)| 2 is drawn 
in Fig. 4(b) where \ip(NAx, 0)/ip(0, 0)| 2 ~ exp(— A). Fig. 4(b) is correct only when the 
terminal point is R = NAx. Fig. 4(b) does not serve for H ^ Hn by a simple shift of the 
terminal point away from NAx. 

Under the condition of Euclidean resonance H = Hr the function \i/j(x, 0)| 2 in Fig. 6 is 
periodic at least with the exponential accuracy provided by the method of trajectories. This 
means that \ip(x, 0)| 2 , besides the periodic oscillations, may have a power law decay. The 
spatial oscillations and not exponentially small tunneling through a long barrier draw the 
analogy with motion in a conventional narrow band periodic structure. Like this motion, in 
our case there is a long distance under-barrier coherence set by the magnetic field. One can 
interpret this as formation of a periodic variable gauge potential. 
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An interference of various paths in the problem of localization in disordered systems 
[22] does not depend on a particular shape of an impurity potential and is determined by 
a mean free path. In contrast to this, the under-barrier coherence strongly depends on a 
potential shape since v(rj) should have a well as in Fig. 2. This rule provides a choice of 
u(y). An arbitrary u(y), for example u(y) = uoy 2 /a 2 , does not result in that well but the 
potential u(y) = u (y 2 /a 2 + y A /a A ) does. Analogously, a pure harmonic potential u{y) = 
uq (1 — cosy /a) is not suitable and should be supplemented by a double harmonic at least. 

Below we consider, as an example, the double harmonic potential u(y) = 
u (1 — cos y/a) (1 — A cosy /a) which is created between two quantum dots on a surface 
with a perpendicular magnetic field. It is sufficient to have this analytical form only at the 
length \y\ < Arj. The dots are separated by the distance R. Sizes of the dots determine 
the distance between their discrete energy levels which is supposed to be \E\ ~ 0.01 eV 
[23]. The parameters of the double harmonic potential are a = 50 A, A = 0.215, and 
u = 1 eV. The spatial period of the potential is 2na = 314 A. After calculations one can 
obtain H R ~ 10 Tesla, Ax ~ 120 A, At] ~ 110 A, and AA ~ 12. At H = H R the tunneling 
probability does not fall exponentially with R which enables to consider tunneling through 
long barriers, say, R = 1 cm. 

The long distance under-barrier coherence can be influenced by dissipation. Dissipation 
results in a finite width 5E of energy levels inside the wells and also disturbs an under- 
barrier motion, according to Caldeira and Leggett [24]. In the classical dynamics dissipation 
corresponds to the form m'x + m r )x. Using the theory [24] , one can obtain (we omit details) 
the criterion O.27T0 < 1 when dissipation does not influence the frictionless motion. Since 
5E ~ frrf this criterion is equivalent to 5E/E < 7AX/R where the parameter A = h/ ^/m\E\ 
can be interpreted as a de Broglie wave length. When SE/E ~ 0.1 [23] one can estimate 
R < 1000 A. For a bigger R dissipation modifies results and this is a matter of a further 
study. A non-homogeneity Su(x) (including applied voltage) of a barrier in the x-direction 
does not violate the above results as soon as 5u(x) < 4\E\. 

In summary, an answer to the question in the beginning of the paper is positive. The 
nature allows a long distance under-barrier motion which is counterintuitive and contrasts to 
WKB-like tunneling. An example of this motion, tunneling between two wells in a magnetic 
field, is considered in the paper. In the absence of dissipation the magnetic field sets a 
long distance coherence under the barrier and a particle can tunnel through a long potential 
barrier as through a conventional narrow band periodic structure. This phenomenon relates 
to Euclidean resonance. 

I thank A. Barone, A. Bezryadin, G. Blatter, M. Gershenson, V. Geshkenbein, L. Ioffc, 
J. Knight, G. Pepe, A. Ustinov, and R. Webb for discussions of related topics. 
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FIGURES 




FIG. 1. A particle passes (the dashed line) through the tunnel potential. 




FIG. 2. The effective potential for the transverse motion forms the well. 



8 




FIG. 3. The classical trajectory in imaginary time. It is chosen N = 3. (a) The transverse 
component (y = irj). (b) The motion in the direction of tunneling. 




FIG. 4. H < Hr; (a) The action A and its WKB part as a function of distance R = NAx, 
shown by dots, between the two wells, (b) Oscillations in the wave function when R = NAx 
(N = 3). The dashed curve shows the WKB-like dependence in the absence of the magnetic field. 
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FIG. 5. The probability of tunneling as a function of magnetic field has the peak at Hr (Eu- 
clidean resonance), correspond to the condition R = NAx. 




FIG. 6. The case of Euclidean resonance H = Hr. The total distance is R = NAx (N = 3). 
The dashed curve represents the WKB-like dependence. 
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